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ABSTRACT 

WLAN localization has become an active research field in 
recent years. Due to the wide WLAN deployment, WLAN 
localization provides ubiquitous coverage and adds to the 
value of the wireless network by providing the location of 
its users without using any additional hardware. However, 
WLAN localization systems usually require constructing a 
radio map, which is a major barrier of WLAN localization 
systems' deployment. The radio map stores information about 
the signal strength from different signal strength streams at 
selected locations in the site of interest. Typical construc- 
tion of a radio map involves measurements and calibrations 
making it a tedious and time-consuming operation. 

In this paper, we present the design, implementation, and 
evaluation of the AROMA system that automatically con- 
structs accurate active and passive radio maps for both device- 
based and device-free WLAN localization systems. AROMA 
has three main goals: high accuracy, low computational re- 
quirements, and minimum user overhead. To achieve high 
accuracy, AROMA uses 3D ray tracing enhanced with the 
uniform theory of diffraction (UTD) to model the electric 
field behavior and the human shadowing effect. AROMA 
also automates a number of routine tasks, such as importing 
building models and automatic sampling of the area of in- 
terest, to reduce the user's overhead. Finally, AROMA uses 
a number of optimization techniques to reduce the computa- 
tional requirements. 

We present our system architecture and describe the de- 
tails of its different components that allow AROMA to achieve 
its goals. We evaluate AROMA in two different testbeds. Our 
experiments show that the predicted signal strength differs 
from the measurements by a maximum average absolute er- 
ror of 3.18 dBm achieving a maximum localization error of 
2.44m for both the device-based and device-free cases. Our 
results also show that ignoring the effect of the UTD in the 
device-free case leads to significant degradation in accuracy 
up to more than 700%. We also relay lessons learned and 
give directions for future work. 
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tion. 

1. INTRODUCTION 

WLANs are installed primarily for providing wireless com- 
munications. However, recent research has shown that WLANs 
can be used in location determination in indoor environments, 
without using any extra hardware 1 3 13. 19 24- 26 1. Acquir- 
ing the location information for a tracked entity unleashes 
the possibility of various context-aware applications includ- 
ing location-aware information retrieval, indoor direction find- 
ing, and intrusion detection. 

There are two classes of WLAN location determination 
systems: device-based, e.g. [3,25] and device-free, e.g. fT3] 
[24) . Device-based systems track the location of a WLAN- 
enabled device, such as a laptop or PDA. On the other hand, 
device-free systems do not require the entity being tracked 
to carry a device and depend on analyzing the effect of the 
tracked entity on the signal strength to estimate the entity's 
position. Device-free localization systems are composed of a 
number of access points (APs) and monitoring points (MPs). 
The MPs, such as standard laptops and other wireless-enabled 
devices, monitor the APs signal strengths and have fixed lo- 
cations. 

Both device-based and device-free systems usually work 
in two phases: an offline training phase and an online loca- 
tion determination phase. During the offline phase, the sys- 
tem collects signal strengths received from different streams 
at different selected locations in the area of interest, and tab- 
ulates them into a so-called radio map. For device-based sys- 
tems, each stream represents the signal strength from an AP 
to the tracked device. For device-free systems, each stream 
represents an (AP, MP) pair and the radio map tabulates the 
effect of the tracked entity on the fixed streams. The differ- 
ence between active and passive radio maps' construction is 
illustrated in Figure [T] 

During the location determination phase, the system uses 
the information stored in the radio map to estimate the user's 
location. Different location determination systems store dif- 
ferent information in the radio map. For example, the Ho- 
rns system (25] stores the signal strength distribution of the 
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Figure 1: Difference between active and passive radio 
maps construction. In a passive radio map, we have a 
histogram per raw data stream, as compared to a his- 
togram per AP. Also, a user does not carry any device 
when constructing the passive radio map. 



signal strength received from each AP, while the Radar sys- 
tem pj stores the average signal strength received from each 
AP. 

Current methods of radio maps' construction use manual 
calibrations making it a tedious and time-consuming opera- 
tion. Furthermore, each time the layout of the environment 
changes or different hardware is used, the whole process of 
constructing the radio map has to be repeated. In addition, 
the process of radio map construction gets more compli- 
cated, in the device-free case, when the number of tracked 
entities increases, since the radio map needs to take all the 
combination of the possible tracked entities' locations into 
account. For example, for a radio map with I locations and 
a system that wants to track up to n entities, the radio map 
needs to store information about Q) possibilities. This em- 
phasizes the need for a method to automatically construct 



the radio maps for an area of interest. 

In this paper we present the AROMA (Automatic genera- 
tion of RadiO MAps) system which can automatically con- 
struct an accurate radio map for a given 3D area of interest. 
AROMA is unique in supporting automatic radio map gen- 
eration for both device-based and device-free localization 
systems. To our knowledge, AROMA is the first system to 
consider radio map generation for device-free systems and 
the first to consider human effects in device-based systems. 
AROMA combines ray tracing with the uniform theory of 
diffraction |14| to model both the RF propagation and hu- 
man shadowing effect. Ray tracing approximates the elec- 
tromagnetic waves as a set of discrete ray tubes that prop- 
agate through the area of interest and that undergo attenua- 
tion, reflection, transmission and diffraction due to the com- 
plexity of an indoor environment. Although ray tracing has 
been used before in site-specific radio propagation predic- 
tion and several tools have been developed, e.g. ||2"} [l"2"||2"0| 
[22][23| , the main focus of such tools was the radio coverage 
problem, i.e. determining the coverage holes given the APs' 
positions. This does not require high accuracy, and there- 
fore, none of these tools account for the human shadowing 
effect on the RF signal. Existing papers on the human shad- 
owing effect, e.g. (7J, illustrate only the theory behind the 
modeling and not its application. On the other hand, propa- 
gation modeling for localization systems requires high ac- 
curacy, where variations in the predicted signal strength can 
lead to large localization errors. Therefore, the AROMA sys- 
tem has three main goals: (1) to automatically construct an 
accurate radio map for a site of interest, (2) to have efficient 
computations, and (3) to incur minimum overhead on the 
user. 

We present the design of the AROMA system and give de- 
tails about its different components and how they interact to 
achieve its goals. We also evaluate the system under two 
testbeds for the device -based and device-free cases. 

The rest of the paper is organized as follows. Section [2] 
presents the details of the AROMA system. We evaluate the 
performance of the system under two different testbeds in 
Section|3] We discuss our experience while building the sys- 
tem in Section |4] In Section we discuss related work. Fi- 
nally, Section[6]concludes the paper and gives directions for 
future work. 

2. THE AROMA SYSTEM 
2.1 Overview 

The AROMA System uses site-specific ray tracing, aug- 
mented with the uniform theory of diffraction (UTD), to pre- 
dict the RF propagation in a 3D site. Figure [2] shows the ar- 
chitecture of the AROMA System. The input to AROMA is 
the 3D model of the site of interest which can be imported 
automatically from CAD tools or drawn using any free 3D 
modeling tool found in the market today such as Google 
SketchUp and Blender. The 3D manipulation of the site 
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Figure 2: AROMA system architecture. 

model is done using the j Monkey Engine (jME) (jTHJ which is 
a high performance open-source Java-based gaming engine. 

Besides the 3D model, the user must provide the site- 
specific configuration. This includes the locations and an- 
tenna characteristics of the APs and MPs. The antenna char- 
acteristics include transmitting power, frequency, maximum 
gain and radiation pattern. The system comes with two pre- 
defined antenna radiation patterns: isotropic and half-wave 
dipole antenna. The user can also define a customized radi- 
ation pattern depending on the hardware used. The user has 
two options to specify the site-specific configuration: 

1 . Insert them manually using the UI tool. 

2. Provide a configuration file with their locations and an- 
tenna characteristics. 

The user has similar options when specifying the locations 
of the radio map cells. 

AROMA comes with a built-in DB of the approximate val- 
ues of the RF propagation properties of common building 
materials such as bricks and concrete. The user has the op- 
tions of using this DB or providing customized values using 
the UI tool. 

After providing the 3D model and the site-specific config- 
uration, the user starts the system. The 3D model is first pre- 
processed to extract the edges in the scene. The Ray Tracing 
Engine is the core of the AROMA system and is composed of 
three modules: the Ray Launcher, the UTD Engine, and the 
Ray Receiver. The Ray Launcher samples the electromag- 
netic waves emitted from the the APs into a set of discrete 
ray tubes covering the area of interest uniformly and each 
having an associated electric field. The ray tubes propagate 
into the environment and undergo reflections, transmissions 
and diffractions. The Ray Tracing Engine handles the inter- 
actions of the ray tubes with the environment. The UTD En- 
gine handles the changes in the electric field associated with 
the ray tubes resulting from these interactions. The contri- 
bution of each tube in the final received signal strength at the 



Algorithm 1 Ray Tracer algorithm 
start preprocessing; 
for all access point AP do 
Queue Q <— launchRaysQ; 
for i = 1 to depth step 1 do 
while ^empty(Q) do 
ray <— dequeue(Q); 
event <— inter action{ray , site); 
process{event); 
utd(event); 

if receivable(receiver, ray) then 

receive(receiver, ray); 
end if 
end while 
end for 
end for 

start postprocessing; 



MPs can be found by the Ray Receiver. The tracing result 
is then processed by the Radio Map Generator to generate 
the radio map. A sampled scene with RF prediction levels 
rendered on its floor can also be generated by processing the 
tracing result by the Rendering Engine. The overall algo- 
rithm used in the AROMA system is briefed in Algorithm [T] 
Each component is explained in details in the next subsec- 
tions. 

2.2 Pre-Processing 

The pre-processor works on the 3D model to extract edges 
in the model and to assign different materials to the faces in 
the model. In addition it constructs data structures that helps 
in speeding up the computations. 

2.2.1 Edge Detection 

The 3D model is loaded in the tool as an array of tri- 
meshes. A tri-mesh is a set of triangles covering possibly 
non-contiguous surfaces. An edge in the model can be spread 
over multiple tri-meshes, which need to be merged for accu- 
rate modeling of wedge diffraction and for efficient compu- 
tations. The edge detection module uses the hysteresis algo- 
rithm [ 1 ] to detect the edges in the 3D model. The algorithm 
consists of two phases. During the classification phase, each 
triangle side is assigned a weight which is the largest angle 
between the normals to any two adjacent triangles sharing 
that side. In the detection phase, a triangle side is considered 
an edge if it passes the hysteresis test, where two thresholds 
are defined. If the weight of the side is greater than the upper 
threshold, then the side is considered an edge. If it is lower 
than the lower threshold, the side is discarded. Other sides 
are considered edges if they neighbors an edge. 

2.2.2 Bounding Capsules 

The pre-processor encapsulates edges and wedges by bound- 
ing capsules (a swept sphere containing the object) and as- 
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signs them unique IDs to efficiently test for ray-edge inter- 
section events. Bounding capsules has an advantage over 
bounding cylinders as a capsule and another object intersect 
if the distance between the capsule's defining segment and 
some feature of the other object is smaller than the capsule's 
radius [21]. Bounding capsules are also used to calculate in- 
tersections with the human bodies in the environment and to 
efficiently diffract rays around them. 

2.2.3 Materials Detection 

To handle different materials efficiently, the pre-processor 
assigns all edges that have the same material the same color. 
This makes the 3D model more efficient for processing by 
the jMonkeyEngine. 

2.3 Ray Tracing Engine 

2.3.1 Ray Launcher 

APs are represented as point sources that emit electromag- 
netic waves having spherical wavefronts. These spherical 
wavefronts are divided into a number of ray tubes that cover 
them entirely and have equal area. Each ray tube is repre- 
sented by a ray located at its center. The rays emitted from 
each AP must experience two forms of uniformity [6|: 

1. Large scale uniformity: to guarantee unbiased cover- 
age of rays in the 3D environment, 

2. Small scale uniformity: to guarantee that the angular 
separation between rays is constant. 

These conditions are satisfied by emitting rays through 
the vertices of an icosahedron whose center is the transmit- 
ter. An icosahedron has 12 vertices and the angular separa- 
tion between rays emitted from its vertices equals 69°. To 
achieve a better angular resolution, the face of the icosahe- 
dron is divided into smaller triangles using a tessellation fre- 
quency N [20]. The rays are then emitted from the vertices 
of the formed triangles and the ray tubes are hexagonal in 
shape as shown in Figure [3] A good approximation for the 
angular separation between rays in this case is [6|: 

69.0° 

a = (1) 



Section 4.4 discusses the values used for the ray tracing 
parameters. 

2.3.2 Ray Tracer 

The complex indoor environment causes the ray tubes to 
change their original direction through either reflection, trans- 
mission, or diffraction. This results in a phenomena known 
as multipath fading where the transmitted signal reaches a 
receiver via multiple paths. 

Each ray is traced for different kinds of interactions with 
the 3D site. A ray incident on an object produces a reflected 
and a transmitted rays as shown in Figure [4] Reflected rays 
satisfy the laws of reflection: 
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Figure 3: Ray launching by tessellating the faces of an 
icosahedron. The transmitter is located at the center of 
the icosahedron. Each icosahedron face is tessellated into 
smaller triangles using a tessellation frequency (N). Rays 
are launched passing through the vertices of the pro- 
duced triangles. 
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Figure 4: Reflection and transmission of rays. A ray in- 
cident on a surface experiences reflection and transmis- 
sion. The transmitted ray suffers also from internal re- 
flection and transmission. 

1 . The incident ray, reflected ray, and the normal to the 
reflecting surface are coplanar. 

2. Angle of incidence equals the angle of reflection. 
And the transmitted rays satisfy the laws of refraction: 

1 . The incident ray, refracted ray, and the normal to the 
reflecting surface are coplanar. 

2. The ratio of the sines of the angles of incidence and 
refraction is equivalent to the opposite ratio of the in- 
dices of refraction (Snell's Law). 
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Figure 5: Wedge diffracion. A ray intersecting with an 
edge produces a set of diffracted rays forming a diffrac- 
tion cone. (3q is the angle between the incident ray and 
the edge, s is the direction of the incident ray, e is the di- 
rection of the edge, and s is the direction of the diffracted 
ray. Q a is the point of intersection, and P is any point in 
the space lying on the produced cone. 

Where 9± is the angle of incidence, 82 is the angle of re- 
fraction, m and ri2 are the refractive indices of the first and 
second mediums respectively. 

A ray incident on an edge produces a set of diffracted rays 
forming what is called a diffraction cone (as shown in Figure 
|5]l obeying the law of diffraction: 

A diffracted ray and the corresponding inci- 
dent ray make equal angles with the edge at the 
point of diffraction, and they lie on opposite sides 
of the plane normal to the edge at the point of 
diffraction. 

sin/?o = \s x e| = \s x e| (3) 

Where f3 is the angle between the incident ray and the 
edge, s is the direction of the incident ray, e is the direction 
of the edge, and s is the direction of the diffracted ray. 

Each time the ray makes an interaction with the environ- 
ment, its depth is incremented. The tracing of a ray ends in 
one of two cases: 

1. The depth reaches a maximum user-defined value. 

2. The power associated with the ray decreases below a 
defined minimum value. 

2.3.3 UTD Engine 

Electromagnetic waves are discretized using Geometrical 
Optics (GO) into a set of ray tubes. Each ray tube propa- 
gates in a direction identified by the ray located at its center, 
and has associated electric and magnetic fields orthogonal 
to each other and to the direction of propagation. GO can 
account for reflection and transmission fl4"| . The reflected 




Reception Sphere Reception Sphere 

Figure 6: Reception sphere. The size of the reception 
sphere should be large enough to intersect only one ray 
from a given wavefront. 

and transmitted electric fields are related to the incident elec- 
tric field using Fresnel's Field Coefficients (41. However, 
GO fails to account for the electric fields in the shadow re- 
gions which occur due to electric field diffraction. UTD ad- 
dressed the GO deficiency which predicts zero electric field 
in human shadowed regions, and thus UTD is used in mod- 
eling wedge diffraction. Another reason for using UTD is 
that it overcomes the drawbacks of the Geometrical Theory 
of Diffraction (GTD), which is an extension of geometrical 
optics that accounts for diffraction ITT). GTD suffers from 
some problems 1 14|, the most serious of them is that GTD 
predicts singular electric fields near the transition regions. In 
UTD, The diffracted electric field is related to the incident 
electric field using UTD diffraction coefficients fl4) . 

The construction of device-free radio maps requires mod- 
eling of human's body effect on RF signals. At microwave 
frequencies and higher, the human body constitutes an im- 
passable reflector for electromagnetic waves. That is, inci- 
dent waves are reflected and diffracted off the body, along 
other interactions with the surrounding environment. Pre- 
vious work in human modeling has shown a strong corre- 
lation between the RF characteristics of the human body 
and a metallic circular cylinder [7] in indoor radio channels. 
Therefore, we use a metallic cylinder to model the human 
body with radius 0.15m, and height 2m J7J. 

All equations related to the UTD Engine are illustrated in 
Appendix [A] 

2.3.4 Ray Receiver 

MPs are represented as point sinks. A ray is received if the 
MP lies within its ray tube. Alternatively, and for efficiency, 
we can consider a ray received if it lies within a certain dis- 
tance from the MP. This distance depends on the size of the 
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Figure 7: An example of the output of the rendering en- 
gine. At the sampling step, an isotropic antenna is vir- 
tually positioned at each sampling point and the overall 
sampling result is then bi-cubically interpolated over the 
whole floor area. 



ray tube at the MP. The reception sphere model is based on 
this observation |20|. A sphere is formed around the MP 



such that at most one ray from a wavefront intersects with 
the sphere, as illustrated in Figure [6] The reception sphere 
radius equals to the radius of the circle circumscribed about 
the hexagonal wavefront of the ray tube. Different rays have 
different reception sphere radii calculated as: 



a I 

7! 



(4) 



Where a is the angular separation between rays and I is 
the total unfolded distance traveled by the ray till the per- 
pendicular projection of the receiver on the ray path. 

2.4 Post-Processing 

2.4.1 Radio Map Generator 

A radio map is constructed by dividing the area of interest 
into a number of locations, each location has a correspond- 
ing cell in the radio map. The user should provide a list of 
the radio map cell locations and the required type of radio 
map, i.e. active or passive radio map (Figure[T|. As another 
option, the user can select a radio map spacing and the sys- 
tem can automatically generate the radio map locations. 

2.4.2 Rendering Engine 

The tracer can sample the site area to predict the signal 
strength across the entire site of interest. The tracing result 
are then color-shaded and layered over the site floor by the 
Rendering Engine. At the sampling step, an isotropic an- 
tenna is virtually positioned at each sampling point and the 
overall sampling result is then bi-cubically interpolated over 
the whole floor area (Figure [7J. 

3. SYSTEM EVALUATION 




Figure 8: Device-based experiment layout. The figure 
highlights the locations of APs and radio map locations. 



Pt 


2mW 


Antennagain(G ma x ) 


Z.OdBi 


Frequency(f) 


2AGHz 


Antenna Type 


Isotropic 



Table 1: APs' configurations for the device-based exper- 
iment. 



Our validation is composed of two experiments: one for 
the device-based radio map generation and the other for the 
device-free radio map generation. We used two Cisco Aironet 
1 130G Series 802.1 1G Access Point and two D-Link DWL- 
G650 NICs. The system is currently implemented in the 
Windows OS. A NIC Query (15) driver that provides an API 
for user-level queries of NDIS fT6| devices is used to col- 
lect the signal strength samples. We collected 60 samples 
for each location in each experiment. Each experiment has 
different configurations that will be illustrated in details in 
the next subsections. 

3.1 Device-based Radio Map Generation 

The experiment was conducted in a typical apartment with 
an area of 700ft 2 . The environment contains furniture and is 
composed of different materials, like bricks, concrete, wood 
and glass. We marked 11 different locations covering the 
entire area as shown in Figure [8] The APs' configurations 
are summarized in Table Q] 

Figure [9] shows the simulated and measured values for the 
two APs. Table [2] summarizes the results. The figure shows 
that the tool gives the same trend as the measurements with a 
maximum average absolute error of less than 3.2 dBm for the 
two APs. This shows that the tool can be used to model more 
complex scenarios of behavior in a device-based setting. 

3.2 Device-free Radio Map Generation 

The experiment was held in the same environment as that 
of the device-based radio map experiment. 44 locations were 
chosen and are illustrated in Figure 10 A new location is 



introduced here, namely location 0, which represents the en- 
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Figure 9: Device-based experiment, measured vs. simu- 
lated RSS. Location ID refers to the IDs in Figure |8j 





API 


AP2 


RMSE 

Average absolute error 
Standard deviation 


4.00 dBm 
3.2 dBm 
3.18 


3.4 dBm 
3.1 dBm 
3.10 



Table 2: Device-based experiment, measured vs. simu- 
lated results. 



Figure 10: Device-free experiment layout. The figure 
highlights the locations of APs, MPs, and radio map lo- 
cations. 





API 


AP2 


MP1 


RMSE 

Avg. Abs. Err. 
Stdev 


1.19 dBm 
0.71 dBm 
1.2 


2.11 dBm 
1.61 dBm 
2.13 


MP2 


RMSE 

Avg. Abs. Err. 
Stdev 


3.02 

2.17 dBm 
3.06 


1.75 

1.21 dBm 
1.77 



Table 3: Device-free experiment, measured vs. simulated 
results. 



Radio map 


Measurements- 
based Classifier 


Simulation- 
based Classi- 
fier 


Device-based 


0.89m 


1.61m 


Device-free 


0.96m 


2.44m 



Table 4: Localization performance (mean distance error) 
for the device-based and device-free experiments. 



vironment without the human. The APs' configurations are 
summarized in Table [TJ 



Figure 11 shows the simulated versus measured RSS for 
the experiment. The results are summarized in Table [3] The 
results show that simulated values are close to the measured 
values with a maximum average absolute error of 2.17 dBm 
for all streams. The figure also shows that the human effect 
is maximum when the person is cutting the line-of-sight, e.g. 
locations 11, 34, and 35 in Figure [Tl(d)| This is captured by 
both the tool and the actual measurements which shows that 
the tool can be used to model more complex scenarios of 
behavior in a device-free setting. 

3.3 Localization Performance 

In this section, we compare the performance of a simple 
nearest-neighbor localization classifier, similar to the Radar 
system [3|, when trained by measured and simulated data. 



We divide the collected samples into two parts with a ratio 
of 2 to 1. The larger part is used as a training data for the 
measurements-based classifier and the other as a testing data 
for both types of classifiers. 

The results are summarized in Table [4] For the device- 
based localization system, the mean distance error calcu- 
lated for the measurement-based and simulation-based clas- 
sifiers are 0.89m and 1.61m respectively. Similarly, for the 
device-free case, the mean distance error calculated for the 
measurement-based and simulation-based classifiers are 0.96m 
and 2.44m respectively. Note that since cross-validation is 
used to evaluate the measurement-based classifier, its accu- 
racy is over estimated. Using an independent test set, the 
accuracy of the measurement-based classifier will be worse, 
while the accuracy of the simulation-based classifier will not 
be affected. Therefore, the difference between the measurement- 
based and simulation-based classifiers will be less. 
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3.4 Effect of Ignoring the UTD 

In this section, we show the effect of not using the UTD 
engine, which is the common practice in previous work in 
the area of RF propagation. Turning off the UTD engine 
disables all the interaction with the human except for the at- 
tenuation effect, i.e. the signal is attenuated by a constant 
amount if the human is obstructing the signal. Figure [12] 
shows the results and Table summarizes them. The results 
show that ignoring the contribution of the UTD engine de- 
grades the performance significantly, up to 706%. This is 
specially true for MP2, where its location makes it signifi- 
cantly affected by the UTD effects (Figure [T0|. 

4. DISCUSSION 

In this section, we discuss different aspects of the AROMA 
system and our experience while developing it. 

4.1 Computational Efficiency 

The AROMA system uses a number of techniques to re- 
duce the computational overhead of ray tracing. Using bound- 
ing capsules, combining edges, and mapping materials to 
colors are all performed at the pre-processing stage to in- 
crease the efficiency of different modules while the ray trac- 
ing engine is running. In addition, the ray launcher's param- 
eters can be tuned for maximum accuracy and performance. 
Finally, the concept of reception sphere is used in the ray 
receiver to further reduce computational requirements. 

4.2 Ease of Use 

Our AROMA system contains a number of modules that 
enhance the user experience and automate common tasks. 
The UI module allows the user to use point-and-click to en- 
ter the locations of APs and MPs. Different wizards are pro- 
vided for the user to enter the APs and MPs parameters, e.g. 
the antenna's pattern, and materials properties. A database 
of different materials' properties, along with a set of prede- 
fined antenna patterns are also included to reduce the over- 
head of entering this information. In addition, the user can 
configure the system to automatically select the radio map 
locations. 

4.3 Prediction Precision 

In a perfect world, the measured RSS should match the 
predicted RSS perfectly. There are different reasons for de- 
viation of the predicted values from the measured ones. Ex- 
tensive work has been done in the characterization of signal 
loss for various materials in indoor environments. However, 
in practice the physical and electrical characteristics of ma- 
terials vary considerably from theory. This variance affects 
greatly the accuracy of signal prediction. In a similar man- 
ner, the human body is made up of water that causes the 
signal to attenuate. The amount of signal loss depends on 
the size of a human body and its orientation that varies con- 
siderably from one person to another. 





API 


AP2 


MP1 


RMSE 

Avg. Abs. Err. 

/" 1 ^ i T f ' 1 / 1 '1 f 1 /" A t"* 

/C LlCgl aUaLlUU 

Stdev 


2.87 dBm 
2.60 dBm 
266% 
1.22 


5.91 dBm 
5.49 dBm 
241% 
2.20 


MP2 


RMSE 

Avg. Abs. Err. 
% degradation 
Stdev 


11.53 dBm 
11.01 dBm 

407% 
3.46 


10.19 dBm 
9.76 dBm 
706% 
2.95 



Table 5: Device-free experiment when the UTD engine is 
not used, measured vs. simulated results. 



In addition, there are always discrepancies between hard- 
ware characteristics and the characteristics included in its 
data sheet. The effect of these discrepancies is inevitable. 
The most influencing discrepancy is that resulting from the 
antennas radiation patterns and orientation. The difference 
between where the radio point location is selected and where 
the human is actually standing is another reason for devia- 
tion. In summary, the reasons of deviation of predicted val- 
ues from the measurements include: 

• Lacking of a perfect 3D model of the environment. 

• Lacking of the precise electrical characteristics of the 
materials in the environment. 

• Hardware discrepancies. 

• UTD assumes diffraction off perfect conductors, yet 
the environment contains diverse dielectrics. 

• Surrounding noise and imprecise human location. 

4.4 Ray Tracing Parameters 

The ray tracing algorithm has two parameters namely: tes- 
sellation frequency of APs (N) and tracing depth (d). Through- 
out our experiments, we noticed that the results are not sen- 
sitive to the tracing depth. The reason behind this is that rays 
at higher depths have smaller power and thus make little con- 
tribution to the received signal. We used a tracing depth of 4 
in our experiment. A tessellation frequency between 4 and 8 
gives both good accuracy and efficient computations. 

4.5 Radio Maps for Probabilistic Localization 

Systems 

The localization results reported in this paper are for de- 
terministic location determination systems that depend on 
storing the average RSS in the radio map. Probabilistic loca- 
tion determination systems, e.g. (25), store RSS distributions 
in the radio map. In order to generate a probability distribu- 
tion, different approaches can be used: Perturbing the loca- 
tion of the device or human, for device-based and device- 
free systems respectively, leads to changing the generated 
radio map and hence can be used to estimate the probabil- 
ity distribution. Similarly, the locations and antenna patterns 
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of APs and MPs and the state of doors and windows can be 
changed to achieve a similar effect. Another possibility is 
to develop models for the relation between the average RSS 
of a stream and its variance. Previous studies have reported 
that the higher the average RSS, the higher the variance (25j. 
This effect can be more accurately quantified and used. 

4.6 Dynamic Radio Map Generation 

The radio map can change also from time to time based 
on the environment conditions. To capture these changes, 
the antennas' parameters and the building materials can be 
adjusted based on environment changes, such as temperature 
and/or humidity over the day or for a given semester in the 
year. Other effects that are dependent on the time of the day, 
such as human density, can also be included in the tool, e.g. 
by randomly placing humans in the area of interest based on 
the time of day. 

Feeding the tool with a small number of tuning measure- 
ments from the environment can be used as another tech- 
nique to change the system's parameters to capture the changes 
in the environment. 

4.7 Human Modeling 

Currently, the AROMA system models the human as a per- 
fectly conducting cylinder |7]. Although this gives good re- 
sults, we believe that the system's accuracy can be enhanced, 
especially for the device-free case, by having a better human 
model. The main drawback of the cylinder model is that it 
does not capture the effect of the user's orientation. This is 
an area for ongoing work. 

4.8 Other Uses 

In this paper, we showed how the AROMA system can be 
used in constructing accurate radio maps for device-based 
and device-free localization systems. However, the AROMA 
system can also be used in a number of different applica- 
tions. In addition to the traditional site planning, such as de- 
termining the placement of APs to maximize coverage and 
detect coverage holes, AROMA can be used to select the lo- 
cations of APs and MPs for maximum combined coverage 
and localization accuracy. AROMA can be also used in the 
design of new algorithms for localization systems, especially 
for the challenging and open research area of device-free 
tracking and identification p4| . The system can be used to 
test the effect of the existence of different humans in the area 
on the constructed radio map and the effect of different in- 
stance of the same object. Another application for the system 
is in studying the effect of the new smart APs that have auto- 
matic power adjustment on the performance of localization 
systems. 

5. RELATED WORK 

Indoor radio propagation has been an active field of re- 
search (2][T2][20 22 23], especially for cellular networks. 
Seidel and Rappaport discuss different radio propagation mod- 



els in |20|. Based on the ray-tracing technique, most au- 
thors model the propagation loss using the simple attenua- 
tion model: 

P(d)=P {d ) - 10 x n x log 10 (5) 

Where P (d) is the power at distance d from the transmit- 
ter, P (do) is the power at a reference distance, and n(> 2) 
is the attenuation exponent. 

The simple attenuation model requires extensive measure- 
ments in order to estimate a best-fit attenuation exponent that 
reduces the overall error for a certain site. In |9), Hills et 
al. report that an attenuation exponent of 2.60 results in the 
best fit in the buildings on the Carnegie Mellon University 
campus. Beside requiring extensive site-specific measure- 
ments to estimate its parameters, the model does not account 
for attenuation effects from typical environments, like doors, 
walls, and other different materials. 

A more sophisticated model is the partition model. The 
partition model takes into consideration path loss caused by 
indoor partitions, like walls and doors. This simple model 
cannot give sufficient accuracy for localization systems. For 
example, the Radar system, which accounts for multiple wall 
attenuation effects along the direct path between transmitters 
and receivers, has a mean error of 3.4m RJ. As we show in 
Section[3]and Table[4] the AROMA system can give accuracy 
of 1.61 m in a similar environment. In addition, the proposed 
tools has a number of other uses as discussed in Section[4] 

Another model that is closely related to the partition model 
is the site-specific model. The site-specific model behaves 
like the partition model, but it needs more site-specific pa- 
rameters such as materials and geometries. Ali et al. (8) in- 
troduced a novel model that relate the average path power to 
the site-specific parameters. However, such systems do not 
consider the human effect and their goal is to study coverage 
and not localization. 

In [ 10 1, Biaz et al. presented ARIADNE, a dynamic device- 
based indoor radio map construction and device-based lo- 
calization system, given a number of actual measurements, 
the system adapts to temporal changes of radio propagation. 
Global attenuation parameters are automatically estimated 
using simulated annealing search. The system requires site- 
specific measurements and does not handle the human effect 
or the device-free case. The results in Section l3~4l show that 
ignoring the human effect can lead to a performance degra- 
dation of more than 700%. In addition, the system does not 
consider the device-free case. 

Different from these systems, the AROMA system is unique 
in supporting automatic radio map generation for both device- 
based and device-free localization systems. To our knowl- 
edge, AROMA is the first system to consider radio map gen- 
eration for device-free systems and the first to consider hu- 
man effects in device-based systems. AROMA follows a 
site-specific approach combined with the Uniform Theory of 
Diffraction (UTD) and multipath fading. In addition AROMA 
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targets high accuracy for the constructed radio maps, which 
is needed for localization system. AROMA's performance, 
in terms of both predicted signal strength and localization 
accuracy, are discussed in Section[3] 

6. CONCLUSION 

This paper introduced AROMA, a system capable of gen- 
erating site-specific radio maps for device-based and device- 
free localization systems. AROMA combines 3D RF prop- 
agation with human body-scattering effects to achieve high 
accuracy. We described the different components of the ARO- 
MOA system and how they interact with each other. 

We evaluated the system in two different testbeds. The 
results show that AROMA can achieve high accuracy with a 
maximum average absolute error of 3.18 dBm. We also eval- 
uated the performance of the system with typical localiza- 
tion systems. Using the radio map generated by the AROMA 
system, we achieved 1.61m mean error for the device-based 
case and 2.44m for the more challenging device-free case. 
Our results also showed that the UTD significantly enhances 
the performance in the device-free case by more than 700%. 

We also relayed our experience gained while developing 
the system. As part of our ongoing work, we are experi- 
menting with different human models. Using the system for 
developing new algorithms for multiple entity tracking and 
identification for the device-free case is a possible research 
area. We are also looking at different techniques for gener- 
ating dynamic radio maps and probabilistic radio maps. 

Our experience with the AROMA system showed that it 
achieves its goals of: 

• High accuracy: through modeling of the human and 
combining geometrical optics with the uniform theory 
of diffraction. 

• Low computational requirements: through the use of 
pre-processing, bounding capsules and reception spheres. 

• Minimal user overhead: through automating common 
tasks, providing default databases for antenna patterns 
and building materials, and providing a friendly user 
interface . 

Moreover, the proposed tool can be used with a wide set of 
different applications including site planning for both cover- 
age and localization, studying new smart APs with dynamic 
power and channel adjustment, and designing new localiza- 
tion algorithms, among others. 
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APPENDIX 

A. RF PROPAGATION 

Our propagation model uses GO to predict the electric fields at different observa- 
tion points. It takes into account the electric field characteristics like phase shift and 
polarization to make the prediction more accurate. The GO is capable of modeling 
direct, reflected and transmitted electric fields only. However, the reflection and trans- 
mission modeled by GO is for perfect electrically conducting (PEC) surfaces only, 
while a typical 3D environment consists mainly of non-PEC surfaces. A derivation of 
an expression for the reflected and transmitted electric fields for a non-PEC surface 
is provided. We assume that all antennas are vertically polarized and that the emit- 
ted waves have spherical wavefronts. Spherical coordinates are used in describing the 
electric field direction and polarization. Most of the equations in this section are taken 
from [14]. 

A.l Electric Field Propagation 
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For a spherical wavefront, both radii of curvature are equal (pi — p 2 — p). The 
expression describing the propagation of the electric field for a general ray tube at a 
distance s is 



EM =13(0) 



— jks 



(6) 



Where 



• E(0) gives the field amplitude, phase and polarization at the reference point 
s = 0. 

• s is the distance along the ray path from the reference point s — 0. 

• e — ^ 3 gives the phase shift along the ray path. 

• k is the wave number and equals ^ 

• p is the principal radius of curvature of the ray tube at the reference point s — 
0. 

The radius of curvature at distance s becomes: 



Ps= P+s (7) 
And the polarization vector remains the same as the reference point. 

A.2 Electric Field Reflection and Transmission 
for non-PEC surfaces 

Since most surfaces in an indoor environment are non-PEC, a fraction of the elec- 
tric field incident on the surface reflects off the surface and the rest transmits through 
it. The ratio of these fractions depends on the characteristics of both mediums as well 
as the angle of incidence of the incident ray. 

Consider a ray tube propagating in the free space from some source. The ray tube 
impacts a smooth surface at a point Q r . This ray tube is completely described by its 
central ray vector s z , principal radii of curvature p\ and p 2 at the selected reference 
point Q s , and its initial electric field is E 1 " (Q s ) at the reference point Q s . 

The incident electric field E l at Q r can be resolved into two components parallel 
(ejj) and perpendicular (ej_) to the plane of incidence. Similarly, the reflected and 
tran smitted electric fields are resolved into two similar components as shown in Figure 

M 

The ratio of the amplitudes of the reflected and the transmitted electric fields is 
given by Fresnel reflection (T) and transmission (T) coefficients. Fresnel coefficients 
depend on the polarization of the incident field. For a field polarized in the direction 
perpendicular to the plane of incidence, Fresnel coefficients are J4J 



E 7 ^ n 2 cos9 1 ' — nicosf?* 

E 1 ^ n 2 cos9 z +nicos0* 

E^ 2n 2 cosO i 

E 1 , n 2 cos9 l +nicos0* 



(X) 



(9) 



While for a field polarized in the direction parallel to the plane of incidence, Fresnel 
coefficients are : 



-^H — JTicos0* +nicos ^ 



2n2Cos l 



(10) 



(in 



E\ nicosf? 1 -\-n 2 cos6 t 
Where 

• ni is the absolute refractive index of the first medium. 

• ?i2 is the absolute refractive index of the second medium. 

• Q % is the angle of incidence. 

• 6 t is the angle of transmission obtained from Snell's Law[2| 

The reflected field can be expressed in terms of the incident field and Fresnel re- 
flection coefficients as 



E r = [E* • e\) r||g]f + \E i ■ e\j T ± e r ± = E\ r„ e\ + E]_r ± e r ± (12) 

Similarly, the transmitted field can be expressed in terms of the incident field and 
Fresnel transmission coefficients as 



E* = (E* • e\) T|, e* + (V • e l ± j T ± e^ = E\T^e\ + E i ± T ± e t ± (13) 
A negative reflection coefficient means an increase in the phase shift by 180° . 



By the summation of the two components using vector sum, the amplitude and the 
polarization of the resultant the reflected or the transmitted field is obtained. 

The incident ray tube is spherical, p\ — p\ — p, and all surfaces are assumed 
planar. The radii of curvature of the reflected and transmitted ray tubes then reduce to 



(14) 



A.3 Wedge Diffraction 



A shadow region is the area that a GO ray can't reach due to the presence of ob- 
stacles in its path. A shadow boundary is a boundary that defines the shadow region. 
|14| shows that there exist shadow regions for direct and reflected rays. GO is not suffi- 
cient to model these regions since it predicts Zero electric field in the shadow regions. 
Our model uses the Uniform Theory of Diffraction (UTD) to model wedge diffraction 
because of its ability to accurately predict the electric field in the shadow regions as 
well as its abi lity to give a non singular solution at the transition regions. 

Figure p4| shows the geometry of a wedge. A wedge consists of two faces, O-face 
and N-face, and an edge. Each face has unit tangent and normal vectors. Tangents 
are directed away from the edge while normals are directed towards the outside of the 
wedge. The inner angle of the wedge is a. 

The plane of incidence is defined by s and e, while the plane of diffraction is 
defined by s and e. 

The incident electric field is resolved into two components parallel (/3 Q ) and per- 
pendicular (jp ) to the plane of incidence. The diffracted electric field is also resolved 
into two components parallel (/3) and perpendicular (tp) to the plane of diffraction. 

A. 3. 1 Diffracted Electric Field 

Using figure|l5| the incident electric field can be expressed as 



e* = e;^ o + e;,0 

and the diffracted electric field can be expressed as 



(15) 



■->(S P0 t (16) 
The diffracted electric field at a distance s on the path of the diffracted ray is given 



by 



E d - D.E\A{s)e- jks 

Where : 

• D is the UTD diffraction coefficient. 

D= \ 
[ -D h 

• A(s) is the spreading factor defined by: 

A(s) = 



s (p + s) 



(17) 



(18) 



(19) 



• p is the radius of curvature of the incident ray tube at the point of intersection 
with the edge. 

A.3. 2 UTD Diffraction Coefficients for PEC wedges 



f K L\L ro 7 L rn :( f>,<p' : /3 ,n^ - D x + D 2 + R 8)h (D 3 + D 4 ) (20) 

Where R s ^h are the reflection coefficients of the wedge at the edge. For a PEC 
surface R s = — 1 and — 1. 

The components of the diffraction coefficients are given by 



-Di = 



7T+ (4>~ <t> ) 



2ll 



F \kL i a + ( 



V)] (21) 



■-(0-0 ) 
2n 



Da 




F \^kL'a~ ((f> — <f> )J (22) 
F [feL r "a + (0 + 0')] (23) 
F \kL ro a~{<f> + 4>')] (24) 
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The functions a± are defined as: 



2n7r/V a 



(25) 



The meaning and values of the parameters included can be obtained from |14| . 

A. 3. 3 UTD Diffraction Coefficients for non-PEC wedges 

A number of heuristics were made to derive an expression for the diffraction co- 
efficients for non-PEC surfaces. One of these heuristics uses Fresnel reflection and 
transmission coefficients as follows |22| 

D s (L\L r °,L r, \<p,<t>\p ,n) = Di + D 2 + r"D 3 + T°D 4 (26) 



X p - 2feL p cos 2 (fl*) 
Where £(< 0) is the Fock parameter. It is defined from |17| : 

5 = -2m (Q r ) cos(9 i ) 



(30) 



(31) 



Where 



L r °,L rn , 



j>',/3 ,n) =D t + D 2 + r^D 3 +T° ± D 4 (27) 



rjj and rjj 1 are Fresnel's reflection coefficients for parallel polarization of the 
o-Face and n-Face respectively. 

and Vj_ are Fresnel's reflection coefficients for perpendicular polarization 
of the o-Face and n-Face respectively. 



B. HUMAN MODELING 

In order to predict the human body- scattering effects in the indoor environment, a 
UTD-based propagation model is used. In the model, the human body is approximated 
with a perfect conducting circular cylinder 1 7 1 of radius 0.20 m, and height of 2 m. At 
microwave frequencies and higher, the human body constitutes a perfectly conducting 
and impassible reflector for electromagnetic waves [7]. 

B.l Geometrical Configuration 

An incident ray on a cylindrical surface may be either reflected or diffracted, ac- 
cording to the location of intersection between the ray and the cylindrical surface. A 
ray that falls tangential to the cylinder grazes the surface in the same plane of inci- 
dence. An extension of the incident ray beyond the point of grazing at the convex 
surface of the cylinder defines the Shadow Boundary (SB), which splits the space out- 
side the surface into the lit and shadow regions |14|. 

Regions in the vicinity of the shadow boundary are called the transition regions. 
It is at the transition regions where the GO and GTD model fails as the incident field 
drops to zero, and the predicted field becomes infinite, and thus UTD was introduced 
|14| . Reflected rays from the cylinder falls in the lit region, while grazing rays detach 
from the surface into diffracted rays and fall in the shadow zone, which is the non- 
illuminated exterior region of the cylinder. A grazing ray path is determined by the 
shortest distance from the ray source to an observation point P s , and creeping on the 
cylindrical surface |14| . 

B.2 Reflection Off a Cylindrical Surface 

The UTD expression of the field reflected off a cylindrical surface from an incident 
electric field is given by |14|: 



E r {p) = E l (Q r )-R sM -A(s)-e- jk3 

Where: 

• E" 1 (Q r ) is the electric field incident to the cylinder at Q, 

• R Sj h 316 tne soft and hard UTD reflection coefficients. 

• A (s) is the 3-D spreading factor of reflection. 



(28) 



-jks 



is the phase shift of the reflected field. 



• s r is the distance between reflection point and observation point. 

T he s oft and hard reflection coefficients are calculated from the following expres- 
sion {l4\ : 



Where m (Q r ) is the curvature parameter at the reflection point. It is calculated 
as[T4): 



^a t (Q r ) j ■ 



(32) 



Where a± (Q r ) is the radius of curvature of the surface in the plane of incidence at 
Q r . Computation-friendly expressions for the transition and the scattering functions 
can be found in |14| , and (5J. 

L p is the 3-D distance parameter, which is calculated from [17 : 



(p\ + s r ) (p* + s r ) 



(33) 



Where: 



• Pi (P2) 1S me ra dius of curvature of the incident wave in (transverse to) the 
plane of incidence. 

• is the radius of curvature of the reflected wave transverse to the plane of 
reflection. 

In our implementation, waves are assumed to be spherical. In that case, p^ — 
P2 — s 1 where s 1 is the total unfolded distance from the transmitter to the point of 
reflection. The reflected wavefront radii of curvature are calculated from |14| : 



1 

Pl,2 



1 

7^2 



(34) 



Where: 



/l,2 = 



cos(6H) 
± 



sin 2 (e 2 ) sin 2 (0i ) 



sin (#2) sin 



4 



(35) 



The terms aiand a 2 are the principal radii of curvature of the surface at the in- 
cidence point. They measure how the surface bends by different amounts in different 
directions at a given point. In the case of a circular cylinder, a± — radius and 
d2 — infinity, hence the previous expression reduces to: 



ai 



1 sin 2 (0 2 ) 
cos 2 ) \ ai 



(36) 



The quantities sin 2 (#1 ) and sin 2 (#2) are calculated from |l4| : 

sin 2 — cos 2 (a) + sin 2 (a) cos 2 (0 l ) 

sin 2 (O2) — sin 2 (a) + cos 2 (a) cos 2 (^ 1 ) 

Where a is the angle between the principal plane of the surface at the point of 
incidence and the plane of incidence as shown in Figure [17] 
The value of a is calculated from |17| : 



(37) 



(38) 



f.! = -s* • Ul 



(39) 



■■-"•■■-""»■ +{«§}] « 



Where: 

• p* (£) and q* (£) are the soft and hard Fock scattering functions. 

• F(Xp) is the UTD reflection transition function which ensures that fields sur- 
rounding the transition regions remains bounded. 

• ( ^ 0) is the transition function parameter, which is defined from |17| : 



-s ■ U 2 



-1 ti 



(40) 



(41) 



Where U\ and U2 are the principal directions of the surface at the point of inci- 
dence. In the case of a cylindrical surface, U2 is parallel to the cylinder generator, and 
U\ is perpendicular to U 2 - 
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Finally, the reflected spherical wave front radii are calculated as: 

Pi + P2 
Pi. 2 = o 



(42) 



When an incident field reflects from a cylindrical surface, the reflected field am- 
plitude varies from that of the incident field. The 3-D spreading factor governs the 
amplitude variation of the reflected field in 3D. It is calculated from |17| : 



A(s) = 



(Pi + s r ) (p£ + s r ) 



(43) 



B.3 Diffraction Off a Cylindrical Surface 

An incident tangential ray at Qi creeps the cylinder surface and detaches at Qi- 
The location of Q2 depends on the location of the observation point so that the total 
distance from the point of incidence to the observation point, including the creeping 
distance, is minimal. 

The field diffracted from a cylindrical surface is given by [17]: 



E d (p) = E l (QO • T Sth -EC-A(s)- e~ ]ks 

Where: 

• E 1, (Qi) is the incident field on the cylinder at Q± . 

• T s h are the soft and hard UTD diffraction coefficient. 



(44) 



A(s) = 



EC = 



is the 3-D spreading factor of diffraction. 



is the conservation of energy term in the case of spherical 



waves, where t is the creep distance along the surface from Q\ to Q2- 
The soft and hard UTD diffraction coefficients are given by |17| : 



,h = -\Jrn(Qi)m(Q 2 ) ■ 



2 

— ■ e 
k 



Where: 

• F{Xd) is the diffraction transition function. 



2m(Q 1 )m(Q 2 ) 



is the transition function argument. 



Ld is the 3-D distance parameter. 



<Qt) j 



is the Fock parameter. 



Transition and Fock functions are computed as described in the previous section. 
The two radii of curvature of the diffracted wave at the observation point are cal- 
culated from |17| : 

Pi = Pb + ^ + l P2 = S ^ ( 46 ) 

Where: 

• t is the creeping distance. 

• s r is the distance between the observation point and the detachment point Q2- 

• p % h is the radius of curvature of the incident wave in the plane transverse to the 
incident plane, which is calculated from ]17|: 



1 

Pb 



sin (a) cos (a) 
P\ P\ 



(47) 



As we assume spherical waves, that is p\ — — s l , the previous equation 
reduces to: 



Pb ' 



■■ Pi 



(48) 



The final radii of curvature of the diffracted ray are then calculated as the average 
of its two radii of curvature. 
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(c) RSS from AP2 by MP1 
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(d) RSS from AP2 by MP2 



Figure 11: Device-free experiment, measured vs. simu- 
lated RSS. Location ID refers to the IDs in Figure [10} 
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(a) RSS from API by MP1 
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(c) RSS from AP2 by MP1 




2 4 6 8 1 1 2 1 4 16 18 20 22 24 26 28 30 32 34 36 38 40 42 44 
Location ID 



(d) RSS from AP2 by MP2 



Figure 12: Device-free experiment when the UTD engine 
is not used, measured vs. simulated RSS. 



Figure 13: Resolution of electric fields. 




Figure 14: Shadow regions. 




Figure 15: Resolution of electric fields. 




Figure 16: Geometrical configuration of the human 
shadowing problem. 
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